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INTRODUCTION 


The results of the so-called energetic approach to fracture with particular attention to the issue of 
energy dissipation due to crack propagation are applied to the case of a crack with cohesive zone. The 
thermodynamic admissibility of subcritical crack growth (SCG) is discussed together with some 
hypotheses that lead to the derivation of SCG laws. A two-phase cohesive zone model for discontinuous 
crack growth is presented and its thermodynamics analyzed, followed by an example of its possible 
application. 
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INTRODUCTION 


Subcritical crack growth (SCG), under both general and cyclic loading, is a phenomenon that has 
been receiving more and more attention during the last forty years. Starling with early investigations 
mainly on fatigue in metals (Refs. 1-9), current research covers a wide variety of materials, especially 
those such as polymers (Refs. 9-13), and ceramics (Ref. 14), that are becoming important in the 
fabrication of composites. From a theoretical standpoint, the problem is that of relating crack growth to 
the load history. In this sense, fundamental understanding has been provided by the energetic approach to 
fracture (Refs. 15-32), that showed (Refs. 15-19) how SCG is strictly related to the rate of dissipation in 
the vicinity of the crack front 


OBJECTIVE: 

TO RELATE CRACK GROWTH 
TO THE LOAD HISTORY 


A CRACK GROWTH LAW 
AND/OR CRITERION IS NEEDED 
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APPROACHES TO THE PROBLEM 


Several theoretical studies in the continuum thermodynamics of fracture have shown that independently 
of the global or local (around the crack tip) constitutive assumptions, a sharp crack with no cohesive zone 
is constrained to evolve according to the Griffith criterion (Ref. 20). Unfortunately, SCGcannot be _ 
described in terms of the Griffith criterion. In the case of SCG, mainly m fatigue, a number of _£row th 
laws are available, although the great majority of them are based on phenomenological observation only. 


• GRIFFITH CRITERION (1920) 


/i0 IF G z G cr 


ORIGINALLY FORMULATED USING AN ENERGY 
BALANCE APPROACH (FIRST LAW) FOR BRITTLE 
SYSTEMS. 


• FATIGUE GROWTH LAWS (SINCE EARLY 
1950’s) 

- CYCLIC LOADING 

- SUB CRITICAL CONDITIONS 
(GRIFFITH CRITERION DOES NOT APPLY) 

- MOST OF THEM ARE ONLY 
PHENOMENOLOGICALLY BASED 
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CURRENT STATE OF RESEARCH 


Modem continuum thermodynamics sees crack propagation like an internal dissipation mechanism. In 
this sense the propagation of fracture can be described by the evolution of a set of convenient kinematic 
state variables, e.g., crack length, whose driving force can be computed directly from the total free energy 
of the body. The application of the thermodynamics with ISV’s is immediate. One of the important 
outcomes of such an approach is the interpretation of a moving crack tip as a moving heat source and the 
subsequent determination of the corresponding near crack tip temperature field. 


• ENERGETIC APPROACH AS A UNIFIED 
APPROACH: 

- FRACTURE STUDIED WITHIN THE 
FRAMEWORK OF CONTINUUM 
THERMODYNAMICS 

- CRACK SURFACE CONSIDERED AN 
INTERNAL STATE V ART ART E : 


G = - 


dijj 

~di 


- CRACK PROPAGATION IS AN 
INTERNAL DISSIPATION MECHANISM. 
IT CAN BE INCLUDED IN 
CONSTITUTIVE THEORIES WITH I.S.V. 

- FORM OF TEMPERATURE 
SINGULARITY AT THE TIP OF A 
RUNNING CRACK 
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IMPORTANT CONTRIBUTIONS 


The present research effort employs many of the results of the modem thermodynamics approach 
fracture. We therefore list some of the most important contributions of this approach. 


THERMODYNAMIC APPROACH 
TO FRACTURE 


• GRIFFITH (1920): CRACK GROWTH CRITERION 
USING THE FIRST LAW OF THERMODYNAMICS 

• CHEREPANOV, G.A. (1967): APPLICATION OF 
CONTINUUM MECHANICS & CONSIDERATIONS 
BASED ON THE SECOND LAW 

• RICE, J.R. (1968): PATH INDEPENDENT 
INTEGRALS IN ELASTICITY; ENERGY RELEASE 
RATE AS CRACK LENGTH CONJUGATE FORCE 

• GURTIN (1979): APPLICATION OF RATIONAL 
THERMODYNAMICS TO A THERMOELASTIC 
SYSTEM WITH A SHARP CRACK 

• NGUYEN (1980-1985): GLOBAL 
THERMODYNAMIC AND DISSIPATION 
ANALYSIS TO FRACTURE 

GENERALIZATION OF THE GRIFFITH CRITERION 

DERIVED BY A DISSIPATION POTENTIAL 

THERMOMECHANICAL SINGULARITY ANALYSIS 



MAJOR PROBLEMS WITH CURRENT METHODS 


The thermodynamic approach to fracture, in the absence of a cohesive zone, derives the Griffith 
criterion as the only possible consequence of the second law. This result is fatigue since fatigue is an 
example of SCG. Another problem in the analysis of cracks with no C.Z. is the loss of weaving of the 
fracture parameter G for almost all material behaviors except the thermoelastic one, thus including special 
behaviors like that of a process zone around a sharp crack. 


• SOME RESEARCHERS HAVE DERIVED 
SUBCRITICAL CRACK PROPAGATION LAWS 
FROM THE FIRST LAW ALONE: 
THERMODYNAMIC ADMISSIBILITY IS 
DISREGARDED. 

• WHEN THE SECOND LAW IS CONSIDERED 
SUBCRITICAL CRACK PROPAGATION HAS 
RF.F.N SHOWN TO BE THERMODYNAMICAL LY 
INADMISSIBLE 

• FOR THE RUNNING CRACK PROBLEM, 
SINGULARITY ANALYSES SHOW THAT G IS 
MEANINGLESS FOR PLASTICITY AND 
VISCOPLASTICITY AND FOR CERTAIN 
VISCOELASTIC MODELS 

• MODELS THAT INCLUDE PROCESS ZONES 
AROUND SHARP CRACKS DO NOT 
NECESSARILY REMOVE THE 
THERMOMECHANICAL SINGULARITY AT THE 
TIP, NOR SOLVE THE ABOVE PROBLEMS. 
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APPROACH USED IN THIS RESEARCH 


The present research effort introduces a cohesive zone into a continuum mechanics model for SCG in 
order to allow for a thermodynamically consistent description of the problem. After postulating the 
presence of a C.Z. ahead of the crack tip, the circumstances under which SCG is thermodynamically 
admissible will be discussed. The assumption leading to the derivation of the traditional form of fatigue 
growth laws is also discussed and a similar form for discontinuous crack growth laws will be obtained. 


• CONTINUUM THERMODYNAMIC 
FRAMEWORK 

- CLASSICAL FIELD THEORY CAN BE 
USED INSTEAD OF NON-LOCAL 
MODELS 

• COHESIVE ZONE 

- ALL THERMOMECHANICAL 
SINGULARITIES ARE REMOVED 

- CRACK TIP HAS A FINITE SIZE 

• SUBCRITICAL CONDITIONS 

- THERMODYNAMICALLY ADMISSIBLE 

- UNIFIED APPROACH TO STUDY 
FATIGUE AND DISCONTINUOUS CRACK 
PROPAGATION. 


61 


BASIC EQUATIONS AND DEFINITIONS 


The analysis prosecuted is basically a global thermodynamic analysis. It consists of deriving global 
thermodynamic statements for the entire structure by interpreting the pointwise governing equations over 
the whole body. 


POINTWISE GOVERNING EQUATIONS 


p«=o.e.-<?. i+ pr 


pj+ 




-p— iO 
T 


o. +p/=0 

Ji J K/, v 


u=u (t kJi T,a’ r ) 

s=s(e ld ,T,a n ) 


(l) 


( 2 ) 


( 3 ) 


( 4 ) 


( 5 ) 
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Note that the constitutive behavior is assumed to be as general as possible through the use of interval 
state variables (at the pointwise level) together with their correspondent solution equations. 


such that 


dh dh 
o. =p ; s =- — 

' dT 


(6) 


where h=h{x,t) is the Helmholtz free energy: 


h = u-Ts 


a n =Q n (e ki ,T,a m ) ; n,m = l,...,N 


( 7 ) 

( 8 ) 


gr-kT. 


( 9 ) 


ALSO LET 


q.r, 

p^rP^ic— 

STRONG FORM OF THE SECOND LAW 

q T 

P 1 ! mir ^0 5 — ^-^*0 

r 1 mic * ^ 

T 2 
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In Figure 1 we have a schematic representation of the system analyzed. The body contains a single 
edge crack which teminates with a cohesive zone characterized by the points a and (3. 



Figure 1 - Crack with a cohesive zone 



Figure 2 - Two-phase cohesive zone model 


64 



DEFINITION OF A CRACK WITH A COHESIVE ZONE 


From a mathematical viewpoint a crack is represented by a line (surface) of discontinuity for the 
various field variables. The cohesive zone is a portion of the crack line (surface) along which a system of 
cohesive forces is acting, and that is also characterized by its own opening displacement 5*. A quantity 
in brackets represents the jump of that quantity across the cohesive zone. At this point it is possible to 
derive global statements for the first and second law and for the dissipation equation. 


QY)=tt(O:0^P(r)} 
c.z. =(x( O : a (0 * C * P (*)} 


o i (U)=o] i y>j=o ji v j 

6 i (c,o^[«<i ; nm,t)=o 


(13) 


GLOBAL FORMS OF THE LAWS OF 
THERMODYNAMICS 


£ 

dt 


(HO 


Jpu dA -J (a ji nu r g i n i ) dS=- J(o J .8.-fa.]v i ) d( 
B S «(0 


(14) 



(15) 


P(0 



mlc T dA-f oiT dA+fqp, dS~ f h,]v, dC 
S s a(t) 


(16) 


WHERE 


S=dB UC 
C*sC-c.z. 


(17) 


65 



DEFINITION OF THE THERMODYNAMIC 
QUANTITIES FOR THE C.Z. 


The cohesive zone is considered a thermodynamic system with its own characteristics. In order to 
discuss such characteristics and write the two laws of thermodynamics for the cohesive zone above, it is 
necessary to define the C.Z. internal energy e, entrophy 9 , temperature fr and free energy \j/. 


jly^const, 


e =■ 


e(CA 


10 , 


a(r)^Pfr) 

f=P(0 


(18) 


! <p O =const, O^C^a(l) 

<P(C,0, «(r)<C<P(t) 

0, C=p(i) 


fr=T*=T- 
VC: tec.z. 


( 20 ) 


9 =e-tph 


( 21 ) 


FIRST LAW FOR THE C.Z. FROM THE GLOBAL 
STATEMENT AND ABOVE DEFINITIONS 

PM PM 

[i dC=( (a i 6 r [?,]v j .) dc (22) 

«M «(f) 


PM PM 

— /« dt+ 2ya = J (afi r [q^v) d C (23) 

«M «M 

LOCAL FORM 

e=o ( 8 r [? j ]v f (24) 
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DISCUSSION ABOUT SECOND LAW AS 
POSTULATED BY GURTIN 


The formulation of the second law for the cohesive zone is not a trivial matter. We therefore adopt the 
statement given by Gurtin in Ref. 26. We then impose a further restriction on the C.Z. behavior so as to 
restrain the C.Z. to act, as a whole, like an actual dissipative system. 


GLOBAL STATEMENT (COHESIVE ZONE 
ALONE) 


P(0 

/<* 


«(<) 



) d(±0 


( 25 ) 


LOCAL STATEMENT 


» 


( 26 ) 


FURTHER RESTRICTIONS 

m 

f [q j]Vj o 

«(0 


THE ABOVE EXPRESSION IS THE DISSIPATION DUE 
TO THE EVOLUTION OF THE COHESIVE ZONE. 
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DISSIPATION ANALYSIS (Being a a Global..) 


In order to properly discuss the dissipation associated with the C.Z. evolution, i.e., crack propagation 
and C.Z. deformation, the thermodynamic force (G-R) work conjugate of the global state variable a must 
be properly characterized. 


BEING a A GLOBAL INDEPENDENT STATE 
VARIABLE, WE CAN WRITE: 


4 >(**, 0 =<{>(**,« ( 0,0 


( 28 ) 


da dt a=cons '- 


( 29 ) 


FIRST LAW: 

m 


«(0 


56 a P<0 

/ —£) d( + (G-R) a = f [q i ]v j d( 

“(0 


( 30 ) 


where 


f P ?5e 

G= f °iir ; R=[ I 1 dC 
o« 8a L da 


( 31 ) 


CRACK ADVANCEMENT RESISTANCE IS A 
FUNCTION OF THE C.Z. THERMODYNAMIC 
STATE. 
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DISSIPATION ANALYSIS (Case 1) 


We will now consider a special type of C.Z. evolution: crack growth with pure translation of the 
cohesive zone. The above assumption is certainly restrictive, but it yields results analogous to that 
obtained in the study of a crack without a cohesive zone. This leads to the following interpretation: a 
running crack with no cohesive zone behaves like a crack with a cohesive zone when the C.Z. is 
constrained to simply translate with the crack tip. 


CASE 1: - PURE TRANSLATION 

- BARENBLATT ASSUMPTIONS 

A=p(r)-<x(r) (32) 


- PURELY ELASTIC COHESIVE 
ZONE 


RESULTS: 


P(o 

(G- 2yJ<x= f [qjv. dCzO (34) 

o(f) 

5(a) 

G= J a. d6 ( (35) 

5(P)=o 


RESULTS ANALOGOUS TO THOSE FOR THE 
CASE OF A CRACK WITHOUT COHESIVE 
ZONE 
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DISSIPATION ANALYSIS (Case 2) 


When a cohesive zone with general behavior, thus with some being of dissipation mechanism, is left to 
evolve without special constraints, we see from die first and second law for the C.Z. that subcritical crack 
propagation, that is a > 0 when G < R, is an admissible phenomenon. 


CASE 2: - ELASTO-PLASTIC COHESIVE ZONE 
- GENERAL DEFORMATION 


iT' 


dilf 

=<p 

dt) 


(36) 


t X =o j df*0 


(37) 


*T- 



(38) 


FIRST LAW BECOMES 

P(0 P(f) 

/ [^]v, dt+(G-R)u = f [qjv. d^O (39) 

«(0 a(t) 


WHERE 


[q^i=o 


38? 

”aT 


. 3<p 96. g c 

0—^=0 — 

dt dt dt 


(40) 


SUBCRITICAL CRACK GROWTH ADMISSIBLE 
WHEN 


P(0 

f dc >(R-G) a (41) 

«(0 

WHEN a = 0 WE HAVE EVOLUTION OF THE C.Z. 
INTERNAL STATE VARIABLES 
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DISSIPATION ANALYSIS (Case 3) 


In general, thermodynamics does not allow to derive evolution laws for the internal state variables, 
thus for the kinematic variables that describe crack propagation and fatigue, some special assumptions can 
be made that allow us to derive a crack evolution law strictly from the first law of thermodynamics. For 
some cases of slow crack propagation, the principle of the minimum entropy production can be evoked, 
thus leading to a certain form of crack growth law. 


CASE 3 - SLOW CRACK GROWTH 

- DISSIPATIVE COHESIVE ZONE 
(ELASTO-PLASTIC) 


ASSUME PRINCIPLE OF MINIMUM ENTROPY 
PRODUCTION HOLDS 

PM 

f dt= 0 < 

«M 


FROM FIRST LAW: 


a 




\ 



R - G 


( 43 ) 


• CYCLIC LOADING 


INTEGRATE OVER A CYCLE 


Aa = 


AA-A Q 


2y ,-G 


M 


( 44 ) 


SIMPLIFIED FORM 


Aa = 


AA 


( 45 ) 
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DISCONTINUOUS CRACK GROWTH 


The analysis presented so far can be easily extended to describe discontinuous crack propagation 
(DCP). With reference to Fig. 2, we present a two-phase cohesive zone model inspired by the 
experimental work by Hertzberg, et al (Ref. 33). Proceeding as in the case of a single phase C.Z. model, 
assuming that the principle of minimum entropy production holds, we obtain an evolution equation for the 
phase separation coordinate ^ that allows to study DCP. 


A VERY SIMPLE 2-PHASE MODEL 
(HERTZBERG, ET AL., 1979) 

P(0 _ P(0 

/ (Vt-™) dO(G„-fi>*(G r R s )5- / [4,1V, dc 
„(,) 5 < * .» 


where 



( 47 ) 


ASSUME a=0 

£>0 ; 0 <G^<R^ 


( 48 ) 


UNDER CYCLIC LOADING 

AA-AQ 


( 49 ) 


SIMPLIFIED 



OTHER DISCONTINUOUS CRACK GROWTH MODELS 

Very few models of DCP have been presented in the open literature. The ones mentioned below are 
those in Refs. 12 and 34. Further study of the thermodynamics of the process is certainly needed. 


• J.G. WILLIAMS, 1977 

• TWO PHASES 

• NOT THERMODYNAMICALLY BASED 


• K. KADOTA & A. CHUDNOVSKY, 1992 

• SINGLE PHASE 

• BASED ON THE PREDICTION OF A 
PROCESS ZONE CRACK RESISTANCE 
DEGRADATION 
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EXAMPLE 


A very simple example of application of the DCP model presented is given. We have assumed that 

G ^h = G ^nax and Hn = X wherc X is the phase transformation energy per unit lengths of transformed 
material. The cohesive zone has been modeled as a two-phase Dugdale zone. 


GEOMETRY: 

- STRAIGHT CRACK IN UNBOUNDED 
MEDIUM 


LOAD: 

- UNIFORM TENSILE STRESS APPLIED 
AT INFINITY 

- CYCLIC, FROM 0 TO T 0 
CASE1: T 0 =0.1MPa 
CASE2: T 0 =0.2MPa 

- PLANE STRESS 

MATERIAL SYSTEM: 

- PS 


PROPERTIES': 

- E=2.2GPa 

- a cz =18.0MPa 

- f=2.0 

- x=30.0J/m 2 


'Botsis et AI, Int. J. Fract., v.33, 1987, pp.263-276. 
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FIGURES 3 AND 4 


In Figs. 3 and 4 we have a schematic of the geometry, load conditions and detailed view of the 
cohesive zone. 



T 


Figure 3 - Example geometry and load conditions 



Figure 4 - Two-phase Dugdale model 




FIGURE 5 


The trend of C.Z. evolution obtained is shown in Fig. 5. It is easy to recognize the discontinuous 
crack growth pattern. 


5.0003E-3 - 1 


Discontinuous Crack Propagation 
Case 1 : a-O.OIm; T 0 -0.1MPa 



Half Crack Length a 

* * * * * C.Z. Phase Separation t 

* 5 5 0 0 C.Z. Tip Coordinate / 6 


4.9999E-3 


10 12 14 


16 18 

Number of Cycles 


20 


Figure 5 - C.Z. Evolution during DCP 
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Half Crack Length (m) 


FIGURE 6 


Figure 6 shows the trend of the growth of a only. 


Discontinuous Crack Propagation 



Figure 6 - Discontinuous crack propagation 



CONCLUSIONS 


In this work a continuum thermodynamic analysis of a crack with a cohesive zone has been presented. 
In particular, the issue of thermodynamic admissibility of subcritical crack growth has been addressed. 
The theory espoused has been applied to the study of DCP and an approximated DCP law has been 
obtained. An example of application of the DCP law is provided. 


• APPLICATION OF THE CONTINUUM 
THERMODYNAMICS APPROACH TO THE 
CASE OF A CRACK WITH A COHESIVE 
ZONE 

• DISCUSSION OF THE ADMISSIBILITY OF 
SUBCRITICAL CRACK GROWTH 

• SUBCRITICAL GROWTH LAWS OBTAINED 
USING THE DISSIPATION EQUATION FOR 
FATIGUE AND DISCONTINUOUS CRACK 
PROPAGATION 
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SIMPLE EXAMPLE OF APPLICATION OF 
THE DISCONTINUOUS CRACK GROWTH 
LAW 



FUTURE WORK 


Further study is certainly necessary, especially toward a better characterization of the constitutive 
equations for the cohesive zone. Also necessary is the coupling of the presented thermodynamic analysis 
with elements of the kinetic theory of fracture in order to obtain more general crack advancement laws. A 
more accurate analysis is needed for the study of DCP together with a stability analysis. Possible 
applications of the theory include the study of delamination in laminated composites, R-toughening in 
ceramics and problems of matrix-fiber interface degradation in fiber reinforced composites. 


• OBTAIN MORE REALISTIC C.Z. 
CONSTITUTIVE EQUATIONS 

• STABILITY ANALYSIS OF THE CRACK 
PROPAGATION EVENT DURING 
DISCONTINUOUS CRACK GROWTH AND 
TRANSITION FROM D.C.P. TO STANDARD 
(SUBCRITICAL) CRACK GROWTH 

• COMPARISON AND COUPLING OF THE 
PRESENT THEORY WITH THE LATEST 
RESULTS OF THE KINETIC THEORY OF 
FRACTURE (SUBCRITICAL CRACK GROWTH 
IN QUASI-PERFECTLY BRITTLE SYSTEMS) 

• POSSIBLE APPLICATIONS: 

• FRACTURE OF POLYMERS 

• DELAMINATION IN LAMINATED 
COMPOSITES 

• FRACTURE OF FIBER-MATRIX 
INTERFACES IN MMC 
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